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We study the Nernst effect due to vortex motion in two-dimensional granular superconductors 
using simulations with Langevin or resistively shunted Josephson-junction dynamics. In particular, 
we show that the geometric frustration of both regular and irregular granular materials can lead to 
thermally driven transport of vortices from colder to hotter regions, resulting in a sign reversal of 
the Nernst signal. We discuss the underlying mechanisms of this anomalous behavior in terms of 
heat transport by mobile vacancies in an otherwise pinned vortex lattice. 
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The Nernst effect - the generation of a transverse volt- 
age when a temperature gradient is applied to a metal or 
superconductor placed in a perpendicular magnetic field 
- has become an important experimental probe of corre- 
lation effects. For example, recent experiments on high- 
T c cuprates 1 and conventional superconducting films 2 
have found a remarkably strong Nernst signal in a wide 
regime above the critical temperature. Being small in 
most ordinary metals, the Nernst effect is naturally at- 
tributed to superconducting fluctuations, cither of Gaus- 
sian nature^— or due to vortex fluctuations ji&Zi In the 
case of granular superconductors, the added complication 
of geometric frustration may significantly affect transport 
properties. Here we show that anomalous sign reversals 
of the Nernst signal can appear in such systems, as the 
magnetic field is varied. 

Let us define a geometry with a perpendicular mag- 
netic field B z and a temperature gradient — V X T induc- 
ing an electric field E y . The Nernst signal and the 
Nernst coefficient v are then defined by 



B z B z (-V X T)' v ; 

In metals the Nernst effect is typically small, being pro- 
portional to particle-hole asymmetry, and v can be ei- 
ther positive or negatived The sign convention adopted 
here conforms with that used in the recent literaturei^— 
In type-II superconductors, the Nernst effect is usu- 
ally much stronger. There, field-induced vortices dif- 
fusing down the applied temperature gradient will gen- 
erate a transverse electric field E = B x v, where the 
drift velocity of the vortices is v = v{— VT), leading to 
E = v\JT x B. The vortex Nernst effect is thus a di- 
agonal response of the vortex current to a temperature 
gradient. Notably, the sign of v is 'positive if vortices 
are driven down the temperature gradient from hotter to 
colder regions. The only way to obtain a negative value 
of v from the vortex motion is then if a situation arises in 
which vortices move from colder to hotter regions, against 
the thermal gradient. A complementary point of view is 
provided by an Onsager relation, relating the Nernst sig- 
nal ejv and the heat current response J® to an applied 
electric current J y , so that Jf* =Tei^J v . It follows that 



J® = Te^GyyB z v Xl which shows that a negative Nernst 
signal (for B z > 0) implies heat flow in the direction 
opposite that of vortex motion. In this Rapid Communi- 
cation we show that such anomalous behavior can indeed 
be realized in granular superconductors. 

Consider first a regular two-dimensional Josephson 
junction array in a magnetic field corresponding to a 
commensurate filling of flux quanta. At low enough tem- 
peratures the vortices will then form a regular lattice 
commensurate with the array. For example, at half fill- 
ing / = 1/2 on a square lattice, the vortices will order in 
a checkerboard pattern. If the density of vortices is low- 
ered slightly below this filling, vacancies are introduced 
into the system, and in absence of any pinning, these will 
be mobile. An applied temperature gradient could then 
produce a drift of these vacancies from hotter to colder, 
resulting in a net vortex flow in the opposite direction 
and consequently a negative Nernst signal. 

We have confirmed this scenario by numerical simula- 
tions in regular arrays (see Figs. Q] - [5] below). Our re- 
sults show that a negative Nernst signal also can appear 
in moderately random Josephson junction networks. We 
have used two different models for the dynamics of the ar- 
rays, (i) Langevin dynamics, and (ii) resistively shunted 
Josephson junction (RSJ) dynamics. The former cor- 
responds to overdamped model- A dynamics,— while the 
latter takes into account current conservation (but ne- 
glects charging effects, i.e., no grain or intergrain capac- 
itance). Previous simulations have been based on time 
dependent Ginzburg-Landau dynamics, 4 which take into 
account fluctuations of the amplitude of the order param- 
eter, and Langevin dynamics^ equivalent to the model 
we use but with different boundary conditions. 

For both Langevin and RSJ dynamics the supercurrent 
flowing between two superconducting grains is given by 

1% = 19, sin^ - 6j - A, ) , Aij — ^ J*' A ■ dr, (2) 

where Ifj is the critical current of the junction, $o = h/2e 
is the superconducting flux quantum, 9i is the supercon- 
ducting phase of grain i, and A is the magnetic vector po- 
tential. We will take A = A cxt + -|^A, where A ext is con- 
stant in time and corresponds to a uniform magnetic field 
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B = Vx A perpendicular to the array, and A = (A x , A y ) 
is time dependent but spatially uniform, describing fluc- 
tuations in the electric field E = —A = — ^A.^ For 
Langevin dynamics the equation of motion is 
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where r]i is a Gaussian white noise with correlations 
(th) = and (m(t)Vj(t')) = (2k B T<y/h)5 ij 8(t - 1'). The 
sum runs over the set A/i of superconducting grains con- 
nected to i. An additional equation describes the dynam- 
ics of the twists A,— 



7aA = 



ij r ji 



c, 



(4) 



0, 



with a time constant 7a = jL x L y and (£(£)} 
(C/.(*)C/(*')> = {2k B T lA /h)5 flv 5{t - t'). Here the sum 
runs over all junctions in the network and Tji = Tj — r^. 

For RSJ dynamics every Josephson junction is shunted 
by a resistor i?, leading to a total current from i to j 
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where 14, is the voltage across the junction, given by the 
ac Josephson relation, 
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The Johnson-Nyquist noise in each resistor obeys 
(m(t)) = and (l$(t)I%V)) = ^{5 ik 6 ft - 
SuSjk)S(t — t'). The equations of motion are obtained 
from current conservation at each grain, and from the 
expression for the average total current 
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This gives a system of coupled differential equations for 
{9i} and A. We assume periodic boundary conditions in 
every direction above, with a fixed average current den- 
sity J cxt . For open boundary conditions the fluctuating 
twist A is redundant and should be set to zero in the 
corresponding direction. 

Temperature enters only via the noise correlations and 
gets a spatial dependence in the presence of a tempera- 
ture gradient. This allows us to calculate the response 
of the system to a temperature gradient. Note that the 
voltage across the system is obtained directly in the sim- 
ulation from E y = —T^-Ay. It is also possible to calculate 
the linear response via a Kubo formula 
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(Ey(t)jQ(0))dt, 
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where the average heat current density is given by 

— x ijAij J Iij, (9) 



Lr.Ln 27T 



(ij) 



with Xji = Xj — Xi and x\^ — (x-i +xj)/2. For Langevin 
dynamics Iij denotes the supercurrent only, while for RSJ 
dynamics it is the total current ([5])iii Since the temper- 
ature is, by necessity, uniform when using the Kubo for- 
mula it is possible to employ periodic boundary condi- 
tions in this case, to eliminate surface effects. Notice 
that the formulation given above is independent of the 
lattice structure. We consider here square, triangular, 
and random lattices. 

One may think of a disordered granular thin film as 
consisting of a random packing of variable sized grains. 
Every grain is connected to each of its neighbors via 
a tunnel junction with a critical current Zy. Thus we 
end up with a randomly connected array of Josephson 
junctions. We model this by generating a random set 
of points with unit density in a square, subject to the 
condition that their separation is larger than some given 
number d m i n . Different values of <i m i n give different levels 
of heterogeneity, with different widths in the distribution 
of grain size diameters. Nearest neighbors are connected 
via a Delaunay triangulation, with the grains as the cor- 
responding Voronoi cells. Some examples are shown in 
Fig.[3]with a heterogeneity varying from 7% to 28%. 

The equations of motion, Eqs. ([3]) and (jU) for Langevin 
dynamics and Eq. ([7]) for RSJ dynamics, are solved nu- 
merically using a forward Euler discretization with a time 
step of At = 0.02 and 0.04, respectively. Note that al- 
though we use a forward Euler scheme to integrate the 
dynamical variables, it is crucial to use a symmetric dis- 
cretization for the heat current J§ ^ For Langevin dy- 
namics the sampling time is set to 20 ■ 10 6 At, after a 
warm-up of 2 • 10 6 Ai, while the corresponding figures are 
10-10 6 Ai and l-10 6 Ai for RSJ dynamics. In addition the 
results are averaged over 64 or more independent runs. 
We consider systems with periodic boundary conditions 
in both directions with sizes up to 160 x 160, but since 
finite size effects are negligible for systems larger than 
20 x 20, only results for this particular system size are 
presented here. The Nernst signal ejv is calculated from 
equilibrium fluctuations using the Kubo formula ([5]). 
while setting J cxt = 0. The validity of Eq. © and the 
discretization used is confirmed by checking that the two 
ways of calculating ejv (Kubo formula and response to 
a thermal gradient) are consistent. We also verify that 
the same result is obtained from the response of the heat 
current to an applied electric current, via an Onsager re- 
lation. Temperature is measured in units of the Joseph- 
son temperature Ej/ks = I c $o/2'xkB, and the Nernst 
signal ejy is given in units of ks/^ej and 27rfcsi?,/$o f° r 
Langevin and RSJ dynamics, respectively. In the major- 
ity of our simulations we use Langevin dynamics, since it 
is computationally less expensive, and gives qualitatively 
the same behavior (see Fig. 2] below) . Unless otherwise 
stated, the results below are for Langevin dynamics. 

Figure Q] shows the Nernst signal ejv as a function of 
filling / = B (A p i) /<I>o for a square lattice ({A p i) is the 
average plaquette area). The different curves correspond 
to different temperatures. At low filling a sharp increase 
culminating in a maximum around / = 0.05 — 0.15, de- 
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FIG. 1: (color online) Nernst signal ejv versus filling fraction 
/ for a 20x20 square lattice at different temperatures T. No- 
tice how the Nernst signal goes clearly negative in the region 
0-4 iS / iS 0.5. Inset: zoom-in of ejv at T = 0.19 around 
/ = 1/3, where the ejv also becomes negative. 

pending on temperature, is observed. This is followed 
by a decrease in the Nernst signal up to half-filling. The 
tilted-hill profile at high temperatures resembles experi- 
mental data of bulk superconductors i 2 - However, for low 
temperatures the curves have significant structure due 
to geometric frustration as the filling is varied through 
different commensurate values. In particular, notice the 
sign reversal of ejy just below half- filling. The inset shows 
a blowup of the region close to another commensurate 
filling / = 1/3, where yet another such a region of neg- 
ative Nernst signal appears, albeit only in a very small 
parameter regime. 

This anomalous sign of the Nernst signal close to, but 
below, commensurate fillings, such as / = 1/3 and 1/2, 
can be connected to the large rigidity (i.e., a relatively 
high melting temperature) of the vortex lattice there. 
This means that as temperature is raised, the vacancy 
defect lattice will melt first, while the vortices remain 
pinned to the underlying lattice. The vacancies can then 
diffuse down the temperature gradient, resulting in an 
opposite net vortex flow. Raising the temperature fur- 
ther will eventually melt also the vortex lattice, restoring 
a positive Nernst signal. It is reasonable to expect that 
also other regions of negative Nernst signal will show up 
in narrow parameter windows at low temperatures, just 
below different commensurate fillings. This scenario of 
melting transitions has been observed in simulations of 
square Josephson junction arrays at / = 5/11 ~ 0.455. 12 
Furthermore, the rich structure of e^r is reminiscent of 
the structure of the resistance as a function of / seen in 
both simulations 13 and in recent experiments^ on square 
Josephson junction arrays. 

For a triangular lattice (see Fig. ^ ejv displays a simi- 
lar behavior as a function of /, but here the structure due 
to geometric frustration effects is even more pronounced. 
The Nernst signal goes once again clearly negative in a 
region just below half-filling, and is strongly suppressed 
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FIG. 2: (color online) Nernst signal ejv versus filling fraction 
/ for a 20x20 triangular lattice at different temperatures T. 
Here ejv shows more structure as a function of /, but again 
becomes clearly negative for / between 0.4 and 0.5. 

around several other fillings, e.g., / = 1/4 and 1/3. The 
relative size of the negative Nernst signal (~ 20%) is more 
or less the same compared to the square lattice case. 

Note that for perfectly regular arrays the Nernst sig- 
nal is periodic, with period one, as a function of fill- 
ing. Furthermore, there is a vortex-vacancy symmetry 
around half-filling, so that e]y(f) — — ejv(l — /), i.e., the 
Nernst signal is naturally negative over large regions of 
0.5 < / < 1 (not shown in the figures). For random net- 
works these properties are absent, and it is not a priori 
clear that the oscillatory behavior with filling persists. 
Figure displays simulation results of the Nernst sig- 
nal versus filling for a couple of different random lattices 
with varying levels of heterogeneity at fixed temperature 
T = 1. As seen, most of the structure found in regu- 
lar lattices is now gone. The same is true also for lower 
temperatures. There is still a sharp increase at low fill- 
ings, reaching a maximum around / = 0.05, followed by 
a smooth decay with increasing /. A negative region ap- 
pears in the most ordered samples (d m i n = 0.8, a = 0.08) 
for fillings 0.4 < / < 1. When increasing the geomet- 
ric disorder by decreasing d min the region gets smaller, 
but it is still visible up to at least d mm — 0.4, a = 0.23. 
As the filling is increased further, a weak oscillatory ten- 
dency can be seen (Fig. HJ). These sign reversals appear 
to be remnants of the natural periodic behavior of reg- 
ular structures, but with an amplitude which is quickly 
damped as filling or disorder is increased. 

For granular superconductors RSJ dynamics should 
give a more realistic microscopic description of fluctu- 
ations compared to the phenomenological Langevin dy- 
namics. In Fig. 2] we compare results obtained using 
Langevin and RSJ dynamics. The curves are essentially 
identical in the interesting parameter regime where frus- 
tration effects are present. The same figure also shows 
the results for a model where the critical currents of the 
junctions are taken proportional to the contact area be- 
tween the grains [or contact length dj:, in two dimensions, 
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FIG. 3: (color online) Top: Nernst signal ejv versus filling 
fraction / at T — 1 for 20 x 20 random lattices with different 
values of the parameter d m i n . Each curve is an average over 16 
disorder realizations. Bottom: Examples of lattice structure 
and size (diameter) distribution of the grains for different val- 
ues of dmin. The grain size standard deviation a is also given 
in each histogram. 

where the dj- 's are the bond lengths of the dual (Voronoi) 
lattice drawn in Fig. [3]. Here the difference is quantita- 
tively larger, but the qualitative features remain. This 
indicates that the geometric frustration dominates the 
Nernst effect and that current conservation and model 
details are less important in this region. 

In conclusion, we have studied the Nernst effect in 
granular superconductors using a phase only model with 
Langcvin and RSJ dynamics. At low magnetic fields 
the Nernst signal displays a characteristic tilted-hill pro- 
file qualitatively similar to experimental findings. «i£ For 



stronger magnetic fields in regular or weakly irregular ar- 
rays, we have found regions of anomalous sign changes of 
the Nernst signal, which translates into sign changes of 
the Nernst coefficient v = ej<sjB z . This is contrary to 
the common belief that the vortex contribution to the 
Nernst coefficient is always positive. A negative Nernst 
coefficient implies a net vortex flow from colder to hot- 
ter regions, and consequently a change in the dominant 
carriers of heat in the system - from vortices to vortex va- 
cancies. Therefore, the Nernst effect offers a unique way 
to probe the nature of heat carriers in superconducting 
structures. We predict that sign reversals of the Nernst 
signal can be seen in experiments on artificial regular 
Josephson junction arrays as well as in granular super- 
conducting thin films at the appropriate magnetic fields. 
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FIG. 4: (color online) Comparison of the Nernst signal ejv 
versus filling fraction / between RSJ and Langevin dynamics, 
and between models with critical current Ifj = I c and Ifj ~ 
djj for a 20x20 random lattice with d m i n = 0.8 at T = 1. 
Each curve is an average over 8 disorder realizations. (The 
RSJ data extends only up to / = 1.) 
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